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Approximate Fitness Functions 



The present invention relates to a method for the evolutionary optimization, to a 
1 0 computer software program for executing such a method as well as to the use of such a 
method for the optimization of aerodynamic structures. 

With reference to figure 1 at first the well-known cycle of an evolutionary algorithm 
will be explained. 

15 

In a step 1 the object parameters to be optimized are encoded in a string called 
individual. A number of such individuals are randomly generated that comprises the 
initial parent generation and the quality (fitness) of each individual in the parent 
generation is evaluated. In a step S2 the parents are reproduced by applying genetic 

20 operators called mutation and recombination. Thus, a new generation is reproduced in 
step S3, which is called the offspring generation. The quality of the offspring 
individuals is evaluated using a fitness function which is the objective of the 
optimization in step S4. Finally, depending on the calculated quality value, step S5 
selects the best offspring individuals (survival of the fittest) which are used as parents 

25 for the next generation cycle if the termination condition in step S6 is not satisfied. 



Before evaluating the quality of each individual, decoding may be needed depending the 
encoding scheme used in the evolutionary algorithm. Note that the steps S2, S3, S4, S5 
and S6 are cyclically repeated until the condition for termination of the optimization in 
30 step S6 is satisfied. 



The algorithm of this evolutionary optimization can be expressed by the following 
pseudo-code: 



Q 



f:\ablaoa\p\pet\24000\24 1 7 8\ep-_d 



Pnnted:25-09-2001 ! 



10 



t:= 0 

encode and initialize P(0) 
decode and evaluate P(0) 
do 

recombine P(t) 
mutate P(t) 
decode P (t) 
evaluate P (t) 
P(t+1) = select P(t) 
encode P (t+1) 
t:=t+l 
until terminate 



15 Such evolutionary algorithms are known to be robust optimizers that are well suitable 
for discontinuous and multi-modal objective functions. Therefore, evolutionary 
algorithms have successfully been applied to mechanical and aerodynamic optimization 
problems, including preliminary turbine design, turbine blade design, multi-disciplinary 
rotor blade design, multi-disciplinary wing platform design and a military airframe 

20 preliminary design. 



25 



Details on the evolutionary algorithms can be found in Dasgupta et al. „ Evolutionary 
Algorithms in Engineering Applications", Springer- Verlag, 1997 and Miettinnen et al., 
^Evolutionary Algorithms in Engineering and Computer Science", John Wiley and 
Sons. 1999. 



30 



One essential difficulty in generation-based approaches to aerodynamic optimization is 
the huge time consumption due to the high complexity of the aerodynamic analysis and 
the large number of evaluations needed in the evolutionary optimization. To address 
this problem, several methods have been developed for constructing approximate 
models. 
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One widely used method in design engineering is the Response Surface Methodology, 
which uses low-order polynomials and the least square estimations. 

A more statistically sound method is the Kriging model, which is also called the Design 
and Analysis of Computer Experiments (DACE) model. In this method, a global 
polynomial approximation is combined with a local Gaussian process and the Maximum 
Likelihood is used for parameter estimation. Artificial neural networks, including 
Multi-layer Perceptrons (MLP) and Radial Basis Function Networks (RBFN) have also 
been employed to build approximate models for design optimization. 



A MLP neural network can be used for building an approximate model. However, it is 
generally difficult to get a model with sufficient approximation accuracy. One of the 
problems is the lack of training data for the neural network, because for many real 
application data collection is a computationally expensive process. This is especially 
15 true when the dimension of the problem space is high. Due to this, the approximate 
model may be of low fidelity and may even introduce false optima. In this case, 
measures need to be taken to guarantee the correct convergence of the optimization 
algorithm when approximate models are used. 



20 Managing approximate models in optimization based on evolutionary algorithms has not 
caught much attention so far. It is known to use a heuristic convergence criterion to 
determine when the approximate model must be updated. The basic idea is that the 
convergence of the search process should be stable and therefore, the change of the best 
solution should not be larger than a user-defined value. An assumption is that the first 

25 sets of data points are at least weakly correlated with the global optimum of the original 
problem, which is not necessarily true for high dimensional systems. 

An approach to coupling approximate models with evolutionary algorithm has been 
proposed in an attempt to balance the concern of optimization with that of design of 
30 experiments. The main idea is to maintain the diversity of the individuals and to select 
those data points that are not redundant for model updating (online learning). In this 
method, the decision of when to carry out the on-line learning of the approximate 
model is simply based on a prescribed generation delay. 
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Most recently it has been proposed to train a neural network with some initial samples 
to approximate the NK model. During evolution, the fittest individual in the current 
population is evaluated on the original fitness function for every 50 generations. This 
5 individual then replaces the one with the lowest fitness in the training set and the neural 
network is retrained. It has been found that the evolutionary algorithm becomes misled 
by the neural network model when the complexity of the original fitness landscape is 
high. The common weakness in the above methods is that neither the convergence 
properties of the evolutionary algorithm with approximate fitness functions (correct 
10 convergence is assumed) nor the issue of model management is addressed. 

In view of the above it is the object of the present invention to reduce the computation 
time and complexity when managing approximate models for fitness functions. At the 
same time the convergence properties of the algorithm should essentially not be 
15 degraded. 

Therefore, a framework for model management with generation based control is 
proposed. Generally, generation based evolution control is more suitable when the 
evaluation of the fitness function is implemented in parallel. The main idea is that the 

20 frequency at which the original computationally expensive function is called and the 
approximate model is- updated should be determined by the local fidelity of the 
approximate model. By local fidelity, the fidelity of the model for the region where 
the current population is located is meant. The lower the model fidelity is, the more 
frequently the original function should be called and the approximate model should be 

25 updated. Since it is hardly possible to build a globally correct approximate model for 
problems with, large dimensionality of the design space, a local model is of more 
practical importance. With this strategy, the computational cost can be reduced as 
much as possible while the correct convergence of the evolutionary algorithm can be 
guaranteed. 

30 

The above object is particularly achieved by means of the features of the independent 
claims. The dependent claims develop further the central idea of the invention. 
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According to the present invention therefore a method for the evolutionary optimization 
with approximate fitness functions is proposed. At first an initial design is set up. The 
initial design is reproduced to create a plurality of offsprings. The quality of the 
5 offsprings is evaluated by means of a fitness function. Finally, the offsprings having 
the highest evaluated quality value are selected for the next generation cycle. 
Selectively, an original fitness function or an approximate fitness function can be used 
for the quality evaluation step. According to the invention the frequency of the use of 
the original fitness function is adaptively adjustable during the optimization process 
1 0 (online adaptation) . 



The approximate fitness function can be refined (updated) on every use of the original 
fitness function. 

15 A covariance matrix of the probability density function adapted during the evolutionary 
optimization can be used for weighting new samples on every update of the approximate 
fitness function. 



The fidelity of the approximate fitness function can be estimated locally on every use of 
20 the original fitness function to adjust the future frequency of use of the original fitness 
function. 

On every use of the original fitness functions the current model error can be calculated 
on the basis of the difference of the result of using the original fitness function and the 
25 approximate fitness function, respectively. The future frequency of use of the original 
fitness function can then be adapted (online) on the basis of the current model error. 

The approximate fitness function can be implemented by means of a neural network 
model. 

30 

According to a further aspect of the present invention a computer software program for 
executing such a method is proposed. 



According to the present invention such a method can be used for structural 
optimization problems, particularly the optimization of aerodynamic bodies, in which 
case the original function can be calculated by a Navier-Stokes solver. 

Further aspects, objects, features and advantages of the present invention will become 
evident for the man skilled in the art when reading the following detailed description of 
an embodiment taken in conjunction with the figures of the enclosed drawings. 

Figure 1 shows a graphical representation of a generation cycle of an 
evolution strategy, 

Figure 2 is a schematic representation of control cycles using an 
approximate model or an original function, 

Figure 3 shows a flow chart for the online adaptation of the use 
frequency of an original fitness function, 

Figure 4 shows an example of the application of the present invention on 
a spline encoded blade construction, 

Figure 5 shows a simulation of blade optimization with approximate 
models, and 

Figure 6 shows an example for the optimization of blades without 
approximate models. 

Figure 2 shows the technique of calculating a fitness evaluation using both approximate 
models (for example implemented by means of neural network models) and the original 
function according to adaptive control cycle. A control cycle is terminated by one 
fitness evaluation using the original function preceded by no, one or a plurality of 
fitness evaluations using the approximate model. Obviously the original function is not 
used for every generation in a control cycle. The control cycle 1 in figure 2 is 
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composed of two generations using the approximate model and is terminated by one 
generation using the original function, whereas the control cycle 2 is composed of one 
generation using the approximate model and is terminated by one generation using the 
original function. In this example at the end of the control cycle 1 it has been decided to 
increase the frequency of use of the original function. 

Figure 3 shows a flow chart of a process of the present invention for online adaptation 
of the use frequency r|(t) of the original function, wherein t is the generation number. 
After the begin the initial design is fixed in step SI. Then an initial control frequency 
r\(0) is fixed in step S7. Depending on the control frequency T|(t) in an evolution 
control step S8 it is decided whether the evaluation of the fitness should be performed 
with the original function (step S9) or with the approximate model (step Sll). Usually 
the evolution control step S8 is designed such that the original function is used in step 
S9 to terminate a control cycle. 

In case the original function is used in step S9 for the evaluation, the data of the result 
of this evaluation are collected in step S10. 
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25 



At the end of a generation cycle (composed of steps S2 to S6 as explained with 
reference to figure 1) it is decided in step S12 whether the current control cycle is 
finished. If not, the processing goes back in a step S14 to the evolution control decision 
step S8. If the end of a control cycle is reached, go to step S13 in which the 
approximate model (the neural network) is updated using the data collected in step S10. 
Based on an estimation of the model fidelity (as will be explained later on) the future 
use frequency of the original function T|(t) is calculated in step S13. Depending on the 
requirements and the result of the evolution steps S9 or Sll, respectively, either the 
final design is output or the processing goes back in step S14 to the evolution control 
decision step S8 to start a new control cycle. 



30 The theoretical background of the present invention will now be explained. 



The canonical Evolution Strategy (ES) operates on an n-dimensional vector x e IR n . 
The main operator in ES is the mutation operator. However, unlike in the canonical 
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genetic algorithm (GA), the mutation operator obtains a direction during search due to 
the self-adaptation of the strategy parameters. The combination of the genotypes from 
different parents, carried out by the recombination operator in ES or the crossover 
operator in GA, usually plays a less significant role or is completely omitted like in this 
5 study. The Evolution Strategy is typically combined with a deterministic selection 
method, either with elitism, i.e., (fi-hX)-ES, or without elitism, i.e., (n,A,)-ES. At the 
same time, the combination with other selection methods, e.g., the EP tournament 
selection is known to be successful. The (|a,A,)-ES is applied here, thus jj. new 
individuals are only selected from the X offspring. The canonical Evolution Strategy 
10 (ES) can be described as follows: 

x(0 = 3c(/-l) + f (1) 

cr, (0 = cr, (* - l)exp(r'z)exp(rz.), V* € {l,...,*} (2) 

15 

x,z<zIR n , z i7 z**N(Q 9 \), z N(0,a(t) 2 ) (3) 



where x is the parameter vector to be optimized and r, r* and <r, are the strategy 
parameters. The values for r and r* are fixed to 

20 




The z,,z and 7 are normally distributed random numbers and a random number 
vector, respectively, which characterize the mutation exercised on the strategy and the 

25 objective parameters. The cr, are also called step-sizes and are subject to 
self-adaptation, as shown in equation (2). Since they define the variances cr, 2 of the 
normal probability distribution, equation (3), for the mutation in equation (1), they 
stochastically define a direction and a step-length of the search process in the 
n-dimensional space. However, arbitrary directions can only be represented by 

30 correlated mutations, i.e., the co variance matrix of the normal distribution should 
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have non-zero off-diagonal elements. Such correlated mutations are realized by the 
Covariance Matrix Adaptation (CMA) algorithm. 

The derandomized Covariance Matrix Adaptation (CMA) differs from the standard 
5 ES mainly in three respects: 

• In order to reduce the stochastic influence on the self-adaptation, only one 
stochastic source is used for the adaptation of both objective and strategy 
parameters. In the derandomized approach the actual step length in the objective 
10 parameter space is used to adapt the strategy parameter. Therefore, the 
self-adaptation of the strategy parameters depends more directly on the local 
topology of the search space. In the case of the adaptation of one strategy 
parameter, this can be written as 



15 



r(0 = a(t - l)«pA(|f| - E[\N(0 9 l\]) 



(5) 



*(/) = x(t - 1) 4- cr(/)z 



(6) 
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z * AT(0,T). 

Therefore, in the derandomized approach the actual step length in the objective 
parameter space is used to adapt the strategy parameter. This results in the 
following, simple, however successful, effect: If the mutation was larger than 

expected (|f |)£[|aT(0,!|]) , then the strategy parameter is increased. This ensures that 

if this larger mutation was successful (i.e. the individual was selected), then in the 
next generation such a larger mutation will again occur, because u(t) was increased. 

The same argumentation holds if (|f |<£[|iV(6J|]) The parameter d is used to regulate 

the adaptation. 



30 • The second aspect is the introduction of the cumulative step size adaptation. 

Whereas the standard evolution strategy extracts the necessary information for the 
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adaptation of the strategy parameters from the population (ensemble approach), the 
cumulative step size -adaptation relies on information collected during ~ successive 
generations (time averaged approach). This leads to a reduction of the necessary 
population size. The cumulation of steps is termed evolution path by Ostermeier and 
Hansen and is formally expressed as: 



The factor c determines the length of the evolution path and the factor c u = *Jc(2 - c) 
is needed for normalization, which can be seen by calculating the variance of s(t) in 
the limit t oo . 

In the CMA algorithm, the full covariance matrix C of the probability density 
function 



is adapted for the mutation of the objective parameter vector. The following 
description of the CMA algorithm follows the one given in N. Hansen and A. 
Ostermeier, „Adapting arbitrary normal mutation distributions in evolutions 
strategies: The covariance matrix adaptation", in Proc. 1996 IEEE Int. Conf. On 
Evolutionary Computation, pages 312-317, IEEE Press, 1996, details of the 
implementation can be found in M Kreutz, B. Sendhoff, and Ch. Igel, „EALib: A 
C + 4- class library for evolutationary algorithms**, Institut fur Neuroinformatik, 
Ruhr-Universitat Bochum, 1.4 edition, March 1999. 

If the matrix B satisfies C = BB T and z ; * N(0 9 \) , then Bz*N(0 9 C). The 
adaptation of the objective vector is given by: 



s(t) = (l~c)s(t-l) + c u z 



(7) 




(8) 



x(0 = x(t - 1) + 8{t - \)B(t - \)z , z% m N(0,l) 



(9) 



-~ — 
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where S is the overall step size. The adaptation of the covariance matrix is 
implemented in two steps using the cumulative step size approach 
(c cov g (O y l)cmdc e (0,1)) determine the influence of the past during cumulative 
5 adaptation). 

s(t) = (1 - c)s(t - 1) + c u B(t - l)z ( 10) 

= 0 - ^av )C(i - 1) + c cov 5(r)5 r (0 (11) 

10 

The next step is to determine the matrix B from C . Since C = BB T is not 
sufficient to derive B 9 the Eigenvectors of C are chosen as column vectors. The 
reason is that in the last step the overall step size 5 has to be adapted. In order to 
perform this adaptation, the expected length of the cumulative vector s s has to be 
15 known. 



Therefore, the variation vector should be N(0 y T) distributed. 

h (0 = (1 - c)s s (t - 1) + c u B s (/ - l)f (12) 
S{t) = 8{t - 1) - exp(4?,(0|| - x a ) (13) 



20 



B s equals B with normalized columns in order for B s (t - l)z to be N(0,l) 
distributed. x n denotes the expectation of the x n distribution, which is the distribution 
25 of the length of an iV(0,l) distributed random vector. 

The CMA algorithm relies on several external parameters which have to be fixed 
manually. 

30 Parallel Implementation of the Evolutionary Algorithm 
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The basic motivation to implement an evolutionary algorithm in parallel is to reduce the 
processing time needed to reach an acceptable solution. This is badly in need when the 
evaluation of the fitness takes a large amount of time. There are several approaches to 
parallelization and the one adopted here is the global parallelization. In this approach 
there is one single population and the evaluation of the individuals is done in parallel. 
For design problems, the evaluation of the individuals usually takes up the 
overwhelming part of the total time consujption, therefore, a sub-ilinear speed-up can 
be achieved if the global parallelization approach is used. The hardware for the 
implementation of parallel evolutionary algorithms can be very different. In our case, a 
network of computers with multi-processors is used. The implementation of the 
parallelization is realized with the Parallel Virtual Machines library. 

EMPIRICAL CONVERGENCE STUDIES AND EVOLUTION CONTROL 

Convergence of the Evolution Strategy with Neural Networks for Fitness Evaluations 

Now the convergence properties of the evolution strategy when a trained multilayer 
perceptron (MLP) neural network is used for fitness evaluations will be evaluated. 

The investigation is carried out on two benchmark problems, the Ackley function and 
the Rosenbrock function. The Ackley function is a continuous, multimodal test 
function, which has the following form: 



where n is the dimension of the function. A neural network is trained to approximate 
the 2D Ackley function. In the simulation, 600 training data are collected from one run 
of evolution with the Ackley function as the objective function. This is to simulate the 
situation that training data are only available from earlier optimization attempts in real 
world applications. The MLP network used has one hidden layer with 20 nodes. 



/, (x) = -20 exp 0.2 J— ^ x f ex P[ ~~ 2 c os27zx t I + 20 + e 
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A false minimum exists near (x u X2) = (l,-4). This can be ascribed to the poor 
distribution of the training samples. Such samples can be used to simulate the situation 
frequently encountered in practice, i.e. sparse and poorly distributed training data. 
5 Recall that for the discussed problems data collection is expensive. If the evolutionary 
algorithm is run with this neural network model, the algorithm converges to the false 
minimum, which is an expected result. 

To show that the identified problem is a general one, simulations are also conducted 
10 oil the Rosenbrock function. 

As) = £Hx?-x M Y+(x,-L0) 2 ] (15) 

1=1 1 1 

For his function, the evolutionary algorithm is expected to easily find a near minimum 
15 but difficult to locate the global minimum. The output of the fitness function spans 
from thousands to zero and therefore it becomes a little more difficult for the neural 
network to learn the mapping. The simulation is conducted on the 2-D Rosenbrock 
function. Similarly 150 samples are generated from one run of evolution on the 2-D 
Rosenbrock function, which is much smoother compared to the Ackley function. 

20 

Comparing it with the true 2-D Rosenbrock function, it is seen that the left ram of the 
2D-Rosenbrock function becomes almost flat due to the poor distribution of the training 
data. No doubt, evolution with this neural network model is vulnerable to serious 
errors, because the global minimum of this approximate model lies in the area where 
25 the true value is very large. Therefore, an evolutionary optimization based on such an 
approximate model will converge to the false minimum. 

Adding Random Samples 

30 One straightforward idea of dealing with the problem is to add some randomly 
generated samples for neural network training, first without considering the cost of data 
generation. In the simulation, 150 random samples are added to the original training 



14 

data, both for the 2-D Ackley function and the 2-D Rosenbrock function. It is found 
that the neural networks are able to learn the main features of the functions and no false 
minima are present. This is confirmed by running the evolution process with the neural 
network model as the fitness function, in which a near-optimum is found for both 
functions. 

Unfortunately, adding random samples does not work for high-dimensional systems. 
To illustrate this, simulations are conducted on the 12-D Ackley function and the 12-D 
Rosenbrock function. In both cases, 450 randomly generated samples are added to 
1000 training data. 450 random samples are added just because in real applications, it is 
impossible to have as much data as needed. Despite that the neural networks have 
achieved a good approximation on the training data, incorrect convergence has occurred 
for both functions. 

Improvement of Convergence with Controlled Evolution 

In the last subsection, it was shown that using additional training samples is not 
effective to deal with the problem of "incorrect" convergence of evolutionary 
algorithms with neural network models. Therefore the concept of evolution control is 
introduced. Two methods are proposed: . 

Controlled individuals: 

In this approach, part of the individuals (t|) in the population (X in total) are chosen and 
evaluated with the real fitness function. If the controlled individuals are chosen 
randomly, it can be called a random strategy. If the best ti individuals are chosen as 
the controlled individuals, it is called a best strategy. 

Controlled generations: 

In this approach, the whole population of r| generations will be evaluated with the real 
fitness function in every X generations, where i\<\. 

Furthermore, when controlled individuals or controlled generations are introduced, new 
training data are available. Therefore, on-line learning of the neural network will be 

14-11-2000 
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applied to improve the approximation of the network model in the region of 
optimization and in turn to improve the convergence of the evolutionary algorithm. 

First the individual based methods for the 12-D Ackley function are investigated. To 
5 determine the number of individuals (r|) that need to be evaluated using the original 
fitness function to guarantee correct convergence, an r| = l,2,...,ll is tested for both the 
random strategy and the best strategy for a (n, X)= (3, 12)-ES. 

All the results are averaged over 10 runs. It can be seen that for the random strategy, 
10 convergence of the evolutionary algorithm is not achieved until ti>9. When 5 out of 12 
individuals are evaluated with the real fitness function, the reported best fitness is close 
to the true fitness. When r| is larger than 7, a near-optimum is found. 

When some of the individuals are evaluated using the real fitness function, new data are 
15 available and on-line learning can be implemented. In the following, it is investigated 
in which way the algorithm can benefit, if these newly available data are used to train 
the neural network on-line. The evolutionary algorithm reports a correct fitness when 
only 2 individuals are evaluated with the real fitness function. When about 50% of the 
whole population are evaluated with the real fitness function, a good near-optimal or 
20 optimal solution is found. 

Similar results have been obtained for the 12-D Rosenbrock function. It is seen that 
when T]>6, the algorithm converges to the correct fitness value and a near-optimal 
solution is found. When on-line learning is introduced, the convergence properties 
25 improve considerably and the resulting solution is near-optimal when r\>5. In both 
cases, the results are based on an average of 10 runs. 

It is shown that the best strategy works much better than the random strategy and 
on-line learning can further improve the convergence properties. 

30 

A FRAMEWORK FOR MANAGING APPROXIMATE MODELS 
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Now a framework for managing approximate models in evolutionary optimization is 
suggested. Taking into account the fact that parallel evolutionary algorithms are often 
used for design optimization, the generation-based evolution control strategy is used in 
this framework. 

5 

Generation-based Evolution Control 

As explained above, generation-based evolution control is one of the approaches to 
guarantee the correct convergence of an evolutionary algorithm when the approximate 

1 0 fitness model has false minima. The idea in the generation based evolution control is 
that in every A. generations there are r\ generations that will be controlled, where rj<X. tj 
should be generally larger than X/2 to ensure the correct convergence when false global 
minima are present in the approximate model. However, an adaptive control frequency 
(y\/X) is proposed so that the calls of the original expensive fitness function can be 

1 5 reduced as much as possible without affecting the correct convergence of the algorithm. 

Determination of Control Frequency 

The higher the fidelity of the approximate model is, the more often the fitness 
20 evaluation can be made using the model and the smaller tj can be. However, it is very 
difficult to estimate the global fidelity of the approximate model. Thus, a local 
estimation of the model fidelity is used. This is feasible because the evolution strategy 
generally proceeds with small steps, i.e., with the normal distribution small mutations 
are most likely. Therefore the current model error is used to estimate the local fidelity 
25 of the approximate model and then to determine the frequency at which the original 
fitness function is used and the approximate model is updated. The frequency is denoted 
by t\I\. When X is fixed, the frequency is solely determined by r\. 

To get a proper r\ 9 heuristic fuzzy rules can be derived as follows: 

30 

If the model error is large, then tj is large. 
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Since the rule system has only one input variable, the input-output mapping of such a 
fuzzy system can be approximately expressed by: 



?](k + 1) = 77 mm + 



E(k) 



ijlmax Vmin) » 



(16) 
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where [x] denotes the largest integer that is smaller than x, tj^ is the maximal rj, 
Tlmax^ and r| min usually equals 1 so that the information on the model fidelity is always 
available. is the allowed maximal model error and E(k) is the current model error 
estimation, k denotes the k-th cycle of X generations. Every X generations is called a 
control cycle. 



15 



The estimation of the current model error is carried out before the next control cycle 
begins. Suppose all the new data in the last ti generations are valid (in aerodynamic 
design, some design may result in unstable fluid dynamics and therefore the data may 
be invalid) and the population size is P, then there will be tiP data in total. Thus, the 
model error is estimated as follows: 



(17) 



20 where y(i) is the true fitness value and y m (i) is the fitness calculated using a feed- 
forward neural network model: 



M ( n \ 



(18) 



25 where H is the number of hidden nodes, n is the number of inputs, w y and v, are the 
weights in the input layer and output layer, and 9(.) is the logic function 



6(z) = 



l + e~' 



(19) 
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The framework for evolutionary optimization with approximate models can be 
summarized as follows (see also the already explained flow-chart of figure 3): 

Begin 

Initialize P(0), ti(0), X, E^, t^, let k=0 
for t=0 to W 

if (t%X)<ri(k).l 

Use the original fitness function 

Collecting new data 

else 

Use the approximate model 

and if 

if (t%X)=XA 

Estimate the local model fidelity E(k) 
Estimate T](k+1), k=k+l 
Select the data for model updating 
Update the approximate model 

end if 

end for 

End 

Recall that the fidelity of the model is estimated locally based on the error information 
from the last cycle. Therefore, X should not be too large. 

Weighted Online Learning Using Covariance Matrix " 

An important issue in on-line learning is how to select new samples for network 
training to improve the model quality as much as possible. A simple way is to use all 
the new samples, or a certain number of the most recent data, if the neural network can 
learn all the data well. Unfortunately, in on-line learning, only a limited number of 
iterations of training are allowed to reduce the computation time. 
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There are a number of methods for data selection in neural network training, which are 
usually called active learning. However, all of these methods rely on sufficient data in 
order to employ methods from statistics. In the problem outlined in this description, 
5 data are sparse and their collection is computationally expensive. At the same time, 
information about the topology of the search space and even more importantly about the 
direction of the search process is contained in the covariance matrix which is adapted 
during the evolutionary process. To exploit this information, the shape of the normal 
distribution, equation (8), is used to weight the newly generated samples. In this way, 
10 the neural network will put more emphasis on the data points that the evolution strategy 
will most probably visit in the next generation. 

Suppose there are N new samples in total, then the cost function for weighted learning 
is given by 
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e = jjZpiOiyV-y^)) 2 (20) 



where, p(i) is the weight for sample i 
20 p(i) = expQx(/) r C- l x(/)) (21) 

The covariance matrix C is the one which is adapted during the evolutionary search in 
the CMA algorithm outlined further above. Before applying the weights to neural 
network learning, they need to be normalized 
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P(0 = Z~ (22) 



where is the maximal weight among p(i), z = l,2,...,M In this way, the neural 
network is able to learn the most important samples, and those with a weight smaller 
30 than a given threshold are discarded in learning. 
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Application of the present invention: Aerodynamic Design Optimization 

One of the main difficulties in aerodynamic design optimization is the tremendous time 
5 consumption in quality evaluation. Generally, a 2-D or 3-D Navier-Stokes solver with 
turbulence model is used for computing the fluid dynamics, which usually takes hours 
of CPU time to get a solution. 

Table 1: Generation-based control with weighted online learning 

10 



No. 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Ave 
rage 


Best 
Fitness 


20.6 


23.9 


22.9 


16.7 


25.2 


18.1 


36.0 


56.8 


38.0 


22.5 


25.8 


Number of 
calls 


996 


948 


936 


900 


960 


960 


876 


1164 


852 


888 


948 



The objective of the optimization is to maximize the efficiency of a turbine blade and to 
minimize the bias of the outflow angle from a prescribed value. The maximization of 
efficiency is realized by the minimization of the pressure loss. In addition, mechanical 

15 constraints must be satisfied concerning stability and manufacturing. In order to 
describe the two dimensional cross section of the airfoil, a spline encoding based on the 
Non-Uniform Rational B-Splines is used. The spline is constructed from a set on N 
four-dimensional control points (x, y, z, w) that define the control polygon, where x, y, 
z are the three dimensional coordinates and w is a weight for the point. A two 

20 dimensional cross section is obtained by setting the z-coordinate to zero. To minimize 
the dimension of the search space the weights of the control points are fixed so that the 
representation used is a regular B-spline in this example. 

Therefore every control point is described only by two parameters, i.e., the x and y 
25 coordinate of the point. Figure 4 illustrates a blade (solid line) that is generated by a 
spline with seven control points. The dotted line shows the corresponding control 
polygon. The population is initialized with a given blade to reduce the computation 
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time. Two neural network models are used to approximate the pressure loss and the 
outflow angle respectively. The development of the pressure loss and the outflow 
angle during the evolution process are shown in Fig 5. In 185 generations, the 
Navier-Stokes-Solver has been called 874 times, which corresponds to about 73 
5 generations of evolution on the original fitness function only. Then the evolution has 
been run for 73 generations and the results are given in Fig. 6. The pressure loss and 
outflow angle are 0.110 and 69.52 with the approximate models and 0.113 and 69.51 
without the approximate models. The final optimization results with approximate 
models are comparable to those without approximate models, but the evolution process 
10 has converged in about 120 generations, in which fewer than 874 calls of Navier- 
Stokes solver are needed. This achieves an over 30 % reduction of computation time. 
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Claims: 



1 . Method for the evolutionary optimization, 
comprising the following steps: 
10 - setting up an initial population as parents (SI), 

- reproducing the parents to create a plurality of offsprings (S2, S3), 

- evaluating the quality of the offsprings by means of a fitness function (S4), 
wherein selectively a original fitness function or an approximate fitness function can be 
used (SS), and 

15 - selecting the offspring(s) having the highest evaluated quality value as parents (S5), 
characterized in that 

the frequency of the use of the original fitness function is adaptable (S13) during the 
optimization process. 

20 2. Method according to anyone of the preceding claims, 
characterized in that 

the frequency of use defines a control cycle terminated by the use of the original 
function, wherein the approximate fitness function is updated (SI 3) at the end of every 
control cycle. 



3. Method according to claim 2, 
characterized in that 

a covariance matrix of the probability density function adapted during the evolutionary 
optimization is used for weighting new samples at the end of every control cycle. 



25 
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4. Method according to anyone of the preceding claims, 
characterized in that 
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the fidelity of the approximate fitness function is estimated locally at the end of a 
control cycle (SI 3) based on every use of the original fitness function (S9) to adjust 
(SI 3) the future frequency of use of the original fitness function. 

5. Method according to anyone of the preceding steps, 
characterized in that 

it comprises the following steps (SI 3) at the end of each control cycle: 

- calculating the current model error on the basis of the difference of the result of using 
the original fitness function and the approximate fitness function, respectively, and 

- adapting the future frequency of use of the original fitness function on the basis of the 
current model error. 

6. Method according to anyone of the preceding claims, 
characterized in that 

the approximate fitness function is implemented by a neural network model. 

7. Method according to anyone of the preceding claims, 
characterized in that 

in case of the use of the original fitness function the quality of all offsprings of a 
generation is evaluated (S5). 

8. Method according to anyone of the preceding claims, 
characterized in that 

in case of the use of the original function the quality of only a selected number of 
offsprings of a generation evaluated. 



9. Computer software program for executing a method according to anyone of the 
preceding claims. 

10. Use of a method according to anyone of claims 1 to 8 for the optimization of spline 
encoded problems. 
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1 1 . Use of a method according to anyone of claims 1 to 8 for the optimization of 
aerodynamic bodies, wherein the original function is calculated by a Navier-Stokes- 
solver. 
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Abstract 

A framework for managing approximate models in generation-based evolution control is 
10 proposed. This framework is well suited for parallel evolutionary optimization that is 
able to guarantee the correct convergence of the evolutionary algorithm and to reduce 
the computation costs as much as possible. Control of the evolution and updating of the 
approximate models are based on the estimated fidelity of the approximate model. The 
frequency at which the original function is called and the approximate model is updated 
15 is determined by the local fidelity of the approximate model. By local fidelity the 
fidelity of the model for the region where the current population is located is 
designated. The lower the model fidelity is, the more frequently the original function 
should be called and the approximate models should be updated. 

20 (Figure 4) 
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Encode and initialize (t=0) 
Decode and evaluate 



Recombine 



Mutate 




Algorithm 
t: = 0; 

encode and initialize P(0) 
decode and evaluate P(0) 

^ ^° recombine P(t) 
mutate P(t) 
decode P(t) 
evaluate P(t) 
P(t+l) = select (P(t) 
encode P(t+l) 
t:=t+l 
until terminate 



(Offspring) 



S3 



Figure 1: 
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